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Energy injected into the Cosmic Microwave Background at redshifts z < 10^ will distort its spec- 
trum permanently. In this paper we discuss the distortion caused by annihilations of relic particles. 
We use the observational bounds on deviations from a Planck spectrum to constrain a combination of 
annihilation cross section, mass, and abundance. For particles with (s-wave) annihilation cross sec- 
tion, (cr \v\) = ao, the bound is / (mx/MeV)"^ (cro/6 x lO'^'^cm^ s"^) {Q.xx h^Y < 0.2, where 
mx is the particle mass, flxx is the fraction of the critical density the particle and its antipar- 
ticle contribute if they survive to the present time, h = Hq/IOO km s~^Mpc~^, Ho is the Hubble 
constant, and / is the fraction of the annihilation energy that interacts electromagnetically. We 
also compute the less stringent limits for p-wave annihilation. We update other bounds on residual 
annihilations and compare them to our CMB bound. 
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I. INTRODUCTION 

Newly proposed particles, especially dark matter can- 
didates, must evade an ever growing array of empirical 
constraints, ranging from bounds obtained in terrestrial 
laboratory experiments (e.g., reference to limits on 
new cooling sources in stars in our galaxy to con- 
straints on the expansion rate of the universe during 
Big Bang Nucleosynthesis (BBN) , when the universe was 
only a few minutes old |^ . In this paper we add another 
constraint to the list, based on the effect of annihilations 
of relic particles on the cosmic microwave background 
(CMB). 

The CMB energy spectrum provides a direct probe of 
the early Universe at redshifts as high as z 2 x 10®. 
Below this redshift, distortions of the spectrum generally 
cannot be thermalized and are observable today. The 
FIRAS (Far Infrared Absolute Spectrophotometer) in- 
strument on the COBE (Cosmic Background Explorer) 
satellite measured the spectrum and found it to have a 
Planck distribution to within a few hundredths of a per- 
cent This measurement places strong upper bounds 
on any energy injection into the CMB after the thermal- 
ization redshift (e.g., references (|J^). 

Decays of relic particles have been considered as a 
source of CMB distortions , but particle annihilations 
have not. Annihilations are typically ignored because the 
classic WIMP (weakly interacting massive particle) dark 
matter candidates have mass mx ^ GeV and their an- 
nihilations "freeze out" at Tp ~ mx/20 > 50 MeV [||, 
long before the time when the CMB becomes vulnerable 
to distortion (at T ~ 0.5 KeV). Freezeout is defined as 
the time when annihilations cease to change significantly 
the number density of the particle; however, annihilations 
continue eternally at some small rate and their products 
can distort the CMB spectrum if they interact electro- 
magnetically. 

The effects of residual annihilations have been con- 
sidered in a few other contexts: Reno and Seckel 
Hagelin, Parker, and Hankcy ||T^, and Fricman, Kolb, 
and Turner [ pT| computed the effect of annihilation prod- 
ucts on the primordial element abundances. Cline and 
Gao and Gao, Stecker, and Cline ||l3| considered 
the possibility of observing directly the photons from 
annihilations at cosmological distances. Bergstrom and 
Snellman [Q, Rudaz and Rudaz and Stecker fl^ ] 
discuss the detectability of a line source from annihila- 
tions to photons within the Milky Way halo. Jungman, 
Kamionkowski, and Griest conclude that halo an- 
nihilations into particles other than photons probably 
cannot be used to place general constraints on particle 
candidates because of astrophysical uncertainties. 

In this paper we compute the energy injected into the 
CMB by annihilating particles as a function of their mass 
and annihilation rate (i.e., the product of cross section 
and abundance squared). We derive constraints on the 
particle properties by comparison with the observed lim- 



its on chemical potential (/x) distortions, and Compton-y 
distortions (§11). We compare these constraints to similar 
constraints obtained from the production of deuterium 
by photodissociation of primordial helium (§111. A), and 
from the diffuse photon background produced after re- 
combination by extragalactic annihilations (§111. B), and 
annihilations in the Milky Way halo (§111. C). 



II. DISTORTIONS OF THE CMB ENERGY 
SPECTRUM 

We consider first the effect of annihilation products 
on the CMB energy spectrum. The distortion of the 
spectrum takes place in two steps: first the high energy 
annihilation products rapidly dissipate their energy into 
the background photons and electrons, and then the low 
energy background evolves more slowly in an effort to 
restore the Planck spectrum. The permanence of distor- 
tions produced after z ~ 10® is simple to understand in 
the following way: A Planck spectrum with a given pho- 
ton number density must have a specific energy density. 
For z ^ 10®, photon nonconserving processes (double 
Compton scattering and bremsstrahlung) are inefficient 
in the background plasma. Therefore, if energy is in- 
jected into the CMB but not the the correct number of 
photons, a Planck spectrum cannot be restored. We now 
discuss in more detail the form of the distortions pro- 
duced in different redshift intervals. 

Down to recombination at z^cc — 1100, photons with 
i?-y ^ 5 KeV quickly cool and produce heated electrons 
by Compton scattering (76 — > 76) or pair production on 
ions (7A^ Ne^e~) (all photons with > 1 KeV can 
cool if z > 3500). The heated electrons quickly dissi- 
pate their energy by inverse Compton scattering on the 
huge number of CMB photons. This process produces a 
distorted spectrum with phase-space distribution 
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where /(x,0) = l/(e^ — 1) is the Planck distribution, 
X = E /T ^ and the Compton-y parameter is assumed to 
satisfy y « \ |jl^. In our case, where there is ample 
time for all of the input energy to be transferred to the 
CMB, the relation between y and the input energy can 
be found by integrating Equation (|l|) to find the energy 
density as a function of y. The result is 4 y = 5p^/ p-^^ 
where 5p-y is the injected energy and p^ is the energy of 
the CMB photons (see reference [19 for a general review 
of CMB distortions or reference |1J] for a more detailed 
discussion of the Compton- ?/ distortion) . Analysis of the 
COBE FIRAS data set gives \y\ < 1.5 x 10"^ §. 

The Compton-y distortion will be preserved to the 
present time if it is produced after zq — 5.4 x lO^uj^^^^ 
(where ojb = Ob/i^/0.02), the redshift above which elas- 
tic Compton scattering would effectively redistribute en- 
ergy between CMB photons, converting a Compton-y 
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distorted spectrum into a Bose-Einstein spectrum with 
distribution 



exp(a:: + /i) — 1 



(2) 



where /i is the chemical potential. Assuming the input 
number of photons is negligible, which will always be true 
in this paper, the chemical potential is fj, — 1.4 5p~f/p^ 
pO| . The FIRAS hmit on this type of distortion is \p\ < 
9 X 10-5 g. 

Equation (||) describes the equilibrium distribution for 
a fixed total number of photons and amount of energy. 
For z > zdc — 2.1 x 10°ujg ' , photons are produced 
effectively by double Compton scattering (ej — » 677) 
so a Planck spectrum = 0) can be restored for ar- 
bitrary energy input ||20|| . Production of photons by 
bremsstrahlung is already ineffective at z ^ zdq (for the 
observed baryon density). 

To summarize: annihilations occurring at z ^ 2.1 x lO'' 
will be unobservablc, annihilations in the range 5.4 x 
10^ ^ z < 2.1 x 10^ will produce a Bose-Einstein spec- 
trum with chemical potential /i, and annihilations in the 
range 1100 ^ z ^ 5.4 x IC* will produce a Compton-y 
distortion. Annihilations at z < Zicc will not significantly 
affect the CMB energy spectrum, but can be observable 
in the diffuse photon background. 

We now compute the fractional energy injection 
Sp^/pj (where p^ is the energy density of the CMB) 
from annihilations of particle species X. We will assume 
throughout this paper that particle X and antiparticle 
X are not identical (we discuss below how our final con- 
straints are strengthened in the case where the particle 
is its own antiparticle). We also assume that nx = nx, 
where nx (nx) is the number density of X (X) (if there 
is a significant asymmetry the relic density of the less nu- 
merous particle will usually be negligibly small |^ , p^ ). 
The energy produced per annihilating particle is 



f mx 



(3) 



where mx is the mass of the particle and / < 1 is the 
mean fraction of the annihilation energy that interacts 
electromagnetically, and the rate of annihilations per unit 
volume is 



{a-a \v\) nxn 



X J 



(4) 



where (ctq \v\) is the cross section. A useful parameter- 
ization for the cross section is (tTa \v\) = (Jq (T/mx)^, 
where n = for s-wave annihilation (e.g., Dirac neutri- 
nos) and n — I for p-wave annihilation (e.g., Majorana 
particles annihilating into much lighter fermions, see ref- 
erence p3|| ). We integrate the energy injection rate from 
time ti to t2 as follows: 
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where Pann = 2 EaTa is the rate of energy injection 
by annihilations, p-yo is the present energy density in 
the CMB, To is the present CMB temperature, nxo = 
nx (z) (1 + z)-'^, and t^, = 2.4 x 10-^^s. We have assumed 
the universe is radiation dominated (this assumption will 
be accurate enough at the redshifts relevant to our cal- 
culation) so that t ~ 0.301 g^^^^mpi T^^ ee U{1 + z)'^, 
where mpi is the Planck mass, — 3.36 is the number 
of effectively massless degrees of freedom at T << MeV, 
and ti, is fixed by the CMB temperature measured at 
z = 0. The result for n = is 



and for n = 1 
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and Ti = T{zi). 

We use zi ~ zdc and Z2 — zq to find the chemical po- 
tential p = 1.4 Sp^l p-y (this assumes a negligible increase 
in the number of photons, which is valid even in the case 
of an electromagnetic cascade), with the results: 



p = 5.1 A = 2.9 x 10"^/ 



mx 



MeV/ 



(for n = 0), (9) 



and 



/i = 2.0 X 10^ A— = 2.7 X 10"*/ 



mx 
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(for n ^ 1), (10) 



where VLxxh^ — 94 (mx/MeV) {nxo + "^xo) /cm~^. 

Observationally < 9 x 10"^ [Q so, for n = 0, we 
have the bound 



VMeV/ 



6 X 10-2''cm-^s-i 



in XX h^)^ 



< 0.3 



(11) 



Similarly, for n = 1 we find 



^ VMeV/ 



6 X lO-^Ws-i /"""^ 

< 3.3 X 10^ 
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We have chosen to scale the cross section as (cro/6 x 
10~^^cni'^s~^) because this is approximately the value 
at which a non-relativistic particle will have relic density 
Uxx^^ — 1 if it freezes out as a result of its annihilations 
while in thermal equilibrium. The reader should keep in 
mind that Six is not necessarily the contribution of X 
to the present energy density if, for example, the particle 
decays subsequent to distorting the CMB. 

Note that the constraint for n = 1 is roughly mx /Ibc 
times weaker than the constraint for n = 0, where Tdc is 
the temperature at zdc- Constraints on n = 1 particles 
that fall out of kinetic equilibrium at Tk > Tbc can be 
roughly estimated from the given n — I constraints by 
multiplying by a factor Tk/Tdq (see reference p^). 

The Compton-y distortion can be obtained by chang- 
ing the limits of integration in Equation (|^) from zdc 
and zc to zq and Zj-cc- For simplicity, we replace the 
lower limit Zroc by the redshift of matter-radiation equal- 
ity, ~ 2.5 X 10^ iloh^ (where SIq is the present density 
in matter), and we ignore deviations from t oc T^^. The 
precise value of this lower limit is not important (because 



of the log dependence) so we assume z, 



For n = 0, 
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3200. 
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0.15m , 



(13) 



while the observational bound is \y\ < 1.5 x 10~^ [Q, 
giving a bound essentially identical to Equation ( pi] ) . For 
n = 1, y ~ 0.005 /i, so the bound will be significantly 
weaker than Equation (p^). 

In the case of n = 0, we predict that the /x and y distor- 
tions will always appear together, with the relationship 
y ~ 0.15 fi. Fixsen et al. Q do not give joint constraints 
on fi and y so we perform our own linear fit to the data 
in their Table 4, using the formula 



where Iq{v) — B^(To) — GQg{iy) is taken from reference 
(along with the necessary error bars), g{i') = i'^B^{9 K) 
is the galactic contamination model used by and 
dSc/dp is the deviation from a blackbody as parame- 
ter p is varied. We constrain jointly the parameters 
AT, Gi, and dp/ p. Based only on statistical errors, 
the 95% confidence upper bound on the energy injec- 
tion is Sp/p < 1.6 x 10~^. Our analysis is complicated by 
the significant systematic errors, 1 x 10~^ and 4 x 10~^, 
quoted by |4| for p and y, respectively. It is not clear 
how these should be combined with the statistical error, 
so we choose arbitrarily to add the two systematic errors 
linearly, and then add the result in quadrature to the sta- 
tistical error, to find a final 95% confidence upper bound 
dp/p < 4.1 X 10~^. We finally obtain the bound 



■' VMeV/ 



6 X 10^-^''cm-^s^i 



(^xx h^)^ 



< 0.2 



(15) 



The combined constraints are plotted as the solid lines 
in Fig. ^ (for rt = 1 we just use the p constraint). 

We truncated Fig. at mx = 1 GeV, where 
hadronic interactions become important; however, the 
CMB bound should continue to apply at higher energies, 
as long as / is computed to account for non-interacting 
annihilation products. Even neutral annihilation prod- 
ucts can contribute to the CMB distortion if they are 
coupled to the plasma in any way, for example, if they 
heat the background protons through hadronic interac- 
tions (see references |^Jl^ for discussions of hadronic in- 
teractions during this epoch). This is not generally true 
for the other constraints that we review in §111. 

These constraints apply to a particle that is distinct 
from its antiparticle; however, to convert to the case 
where particle and antiparticle are equivalent, it is only 
necessary to make the substitutions i^xx ^ ^x, and 
(To/6 X IQ-^'^cm^s-i cro/3 x 10~27(,j^3g-i rpj^g change 
in cross section normalization cancels the increase in the 
annihilation rate at fixed total contribution to the criti- 
cal density (i.e., for the inequivalent case Ta oc nxnx oc 
ri^-^/4, but for the equivalent case Fa oc n\/2 oc n\-/2). 
The scale 3 x 10~^^cm'^s~^ for the cross section is natural 
also because it gives ftxh^ — 1. The same substitutions 
can be used to convert any of the constraints in §111. 



III. COMPARISON WITH OTHER 
CONSTRAINTS ON RESIDUAL 
ANNIHILATIONS 

In this section we compute bounds, in a form similar 
to Equation (pi]), from the production of deuterium by 
photodissociation of "^He (§111. A), from the diffuse pho- 
ton background produced by extragalactic annihilations 
at z < Zrec (§III-B), and from the diffuse photon back- 
ground produced by annihilations in our galaxy (§111. C). 
We restrict our attention to n = because the n = 1 
bound is very weak in all cases. 



A. Photodissociation of He 

At roughly the same time that the CMB energy spec- 
trum becomes vulnerable to distortion by energy injec- 
tion, primordially produced *He nuclei become vulner- 
able to photodissociation by high energy annihilation 
products. At earlier times the nuclei were protected from 
destruction because photons with high enough energy 
{Ej > 20 MeV) to destroy nuclei instead pair produce 
or elastic scatter on the CMB photons (77 e+e^, or 
77 77). This shielding is effective for > rn^/AA T, 
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where the numerical factor accounts for the fact that pho- 
tons in the high energy tail of the CMB spectrum are still 
very numerous (see reference p5[ , and references therein, 
for a full discussion). Once rn^/AA T > 20 MeV, anni- 
hilation products can dissociate ^He , either directly, or 
indirectly through a cascade when mx ^ TOg/44 T. 

Bounds on energy injection can be derived by con- 
sidering the ^He or the D produced when "^He nuclei 
are dissociated, and requiring that the amount created 
is not greater than the observationally inferred primor- 
dial abundances, (bounds from direct photodissociation 
of deuterium extend to slightly lower annihilation en- 
ergy but are significantly weaker [^). Previous analy- 
ses used ^He-I-D because the upper limit on the 
primordial abundance of D was poorly known; however, 
there is considerable uncertainty in the post-big bang 
production and destruction of ^He so we will use D be- 
cause its abundance is now more robustly measured in 
QSO absorption systems j^, with a very conservative 
upper bound no/riH ^ D/H < 3 x IQ-'^ i^^"^^^) [0' 
or, more probably, D/H < 4 x 10^^ (low D) 

Protheroe, Stanev, and Berezinsky presented a de- 
tailed computation of the (redshift dependent) quantity 
of D produced for a given amount of injected energy, 
Nd{z). With this input, the ratio D/H is given by 



D 
H ~ 

For n = this is 



Npjz) 2Ea{aa\v\)nj, ^^ 
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where we have used the results in [ p9[ and evaluated the 
integral numerically. The constraint on annihilations of 
particles with mx ^ 26 MeV (the energy needed to dis- 
sociate '*He into D) is 



-1 

B 



mx 
MeV 



6 X 10~^''cm3s~^ 

< 1.3 (high D), or 
< 0.2 (low D) , (18) 



weaker than the CMB bound if the high deuterium abun- 
dance is used, but similar for the low deuterium abun- 
dance (recall ojb = Ob/i^/0.02 ~ 1). These constraints 
are plotted as the long-dashed lines in Fig. |l|. 

This calculation is not exactly correct unless mx » 
26 MeV because our use of the results in |^ assumes 
that the annihilation products are energetic enough to 
produce a cascade. A more careful calculation would 
weaken the bound slightly for mx ~ 26 MeV pl . 

Our constraint appears to be several orders of mag- 
nitude stronger than the constraint described by Equa- 
tion (9) in reference |26|. They used the observational 



bound {^He + D)/H < 1.1 x 10~^, which, considering 
that ■^He is produced ~ 24 times more efficiently than 
D 1 29 1 , should give a bound somewhat stronger than our 
Equation (|l8|). We believe that the discrepancy is the 
result of a numerical error in their calculation. Our re- 
sult for the {^He + D)/H calculation agrees with 
however, we do not use this result for our annihilation 
constraint because the deuterium constraint should be 
more reliable 



B. Extragalactic 7 Background 

Recently, Kribs and Rothstein |^ used the observed 
7-ray background to constrain late decaying particles. 
Gao, Stecker, and Cline computed the expected 7-ray 
flux from annihilations of a possible lightest supersym- 
metric particle. Here we generalize this annihilation cal- 
culation to match the form of the constraint derived from 
the CMB. Unlike the previous constraints, the constraint 
we compute now only applies to annihilations directly 
to two photons, which typically will be sub-dominant 
for particles with mass greater than the electron mass 
(assuming particle X does not couple directly to pho- 
tons). Therefore, we will only calculate the constraint 
for mx < 0.5 MeV (constraints for general annihilation 
products would be strongly model dependent). 

Temporarily ignoring the possibility that photons pro- 
duced after recombination (at Zj-oc — 1100) are absorbed 
on their way to the observer, the observed energy flux at 
present from extragalactic annihilations to photons is 



dJ 



2 c 



dEdn in 



dt f^mx ^aiz) 



xS[Eo{l + z)-mx] {l + zY 



(19) 



where f-y is the fraction of annihilations that produce 
a pair of photons with E^ = mx- Using dt/da — 
3 toa^/^ /2 = Hq^g^^^ for an Einstein-de Sitter universe, 
and 6 (EoOT^ — mx) — a^S (a — Eo/mx) E^^, we find 



dJ 
dEdfl 



2 c 2 fmx^'^' 



3.9 X 10" 



cm^ s sr 



"^nMevJ [Eo 



3/2 



6 X 10 



i^xxh') 



(20) 



where we have taken n = 0. The photons we are consid- 
ering, with energy E^ < 0.5 MeV, will in fact lose most of 
their energy by Compton scattering if they are produced 
at z > 200 (see reference |3l]), so the maximum redshift 
factor is mx/Eo = 200. 

We derive a bound on annihilations by requiring that 
the predicted photon background is not greater than the 
observed one. After considering observations of the pho- 
ton background at all relevant energies (see references 
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|g2|-g4[), we find that the ASCA data in the energy 
range 0.8 KeV ^ Eo ^ 2.5 KeV provides the best con- 
straint on the annihilating particles that we consider. 
The background in this range is conservatively bounded 



by dJ/dEdn < 0.36 (So/MeV) cm-2s-isr-i Q. 

For 160 KeV < nix < 500 KeV, the best constraint 
is derived from photons produced at z ~ 200. Setting 
Eo = mx /200, we obtain the bound 



■'''{MeVJ 



-1.42 



6 X 10-27c 



< 7.6 



(21) 



The best bound for lower annihilation energies comes 
from observations at Eo — 0.8 KeV: 



(mcv) 



-1/2 



6 X 10^2''cm3s~i 



< 1.3 . (22) 



Figure shows these constraints as the dashed line. 



C. Annihilations in the Milky Way Halo 

The observability of annihilations to photons in our 
own galaxy has been discussed in many papers, including 
p^-p^ . In this subsection we combine the latest observa- 
tional limits on the photon background with the calcula- 
tion by Kamionkowski |35| and Jungman, Kamionkowski, 
and Griest of the flux expected from halo anni- 
hilations to derive a bound in the form of Equation 
(prj). As discussed in the previous subsection, we only 
consider annihilations to two photons of particles with 
mx < 0.5 MeV. 

Jungman, Kamionkowski, and Griest [pT| assume the 
model 



i?2 



(23) 



where liip) 1 depends on the observation angle ip, 
^ ^xX is the contribution to the critical density 
from all kinds of dark matter, and p^j^ ^ 1 is the den- 
sity of dark matter near the solar radius, in units of 
0.4 GeV cm~^. We have assumed px{r) oc pd{t), and 
used the detector energy resolution, ^EjE ~ 0.2, appro- 
priate for the energy bins in Fig. 10 of reference 

The relevant energy range for the halo annihila- 
tions we are considering, 1 — 500 KeV, corresponds 
to a "bump" in the observed spectrum (see p2t). 
To be conservative, we construct a simple bound by 
comparison with the single power law dJ/dEdfl ^ 
0.022 (£'o/MeV)"^-^cm-2 s~^ sr'^ which is an upper 
bound on the energy flux in the full range 1 KeV ^ Eo 
100 GeV Since annihilations to two photons pro- 

duce a line source at Eo = mx, we have the bound 



CTO 



XX' 



< 0.008 , (25) 



which we show in Fig. |l| (the observationally favored 
value for the total mass density is JId/i^ ~ 0.15 p7| , ^ , 
but we use 0.3 as a more conservative upper bound). For 
particles that survive to the present, this bound is gen- 
erally stronger than the bound from distortions of the 
CMB. Note that our calculation is a somewhat rough es- 
timate because the observational bounds are considerably 
lower at some energies than the power law we adopted 
^2| , the energy bin width AE/E is approximate, and we 
used I{ip) = 1 when its value can be higher for observa- 
tion angles near the galactic center [p7| . 

For particles with mx > 0.5 MeV we might consider 
the observation of annihilation products other than pho- 
tons (e.g., positrons); however, reference concludes 
that these signals cannot be used to conclusively rule out 
dark matter candidates because of astrophysical uncer- 
tainties. 



for the dark matter density distribution in the Galaxy, 
where R is the distance of the Sun from the galactic cen- 
ter, a is the core radius, and r is the distance from the 
center of the galaxy. Note that simulations of cold dark 
matter models predict central cusps instead of a core [ p6[ , 
which could lead to enhanced annihilation signals toward 
the center of the Galaxy; however, we want to be conser- 
vative so we do not assume a cusp. Re-writing Equation 
(10.1) of reference |jl^ to include the possibility that the 
particle X does not make up all the dark matter, we flnd 
that the energy flux at Eo = mx is 



dJ 



3.0 



mx 



dEdn •'''cm^ s sr VMeV 



„0.4 

Pd 



m 



6 X 10~2'^cm3s" 



inxxh"? 



(24) 



IV. DISCUSSION 

Some of the parameter space that can be constrained 
by residual annihilations is covered already by other 
kinds of constraints. To put the annihilation constraints 
in perspective we review the primary astrophysical ones. 

BBN gives limits on the expansion rate of the universe 
at T '--^ 1 MeV (sometimes described as a limit on the 
effective number of light neutrinos, e.g., reference [^). 
The expansion rate would be affected by an additional 
particle with mass in the range where our CMB bound 
is most constraining {mx ^ McV), if the particle's num- 
ber density at the time of BBN is equal to the thermal 
equilibrium value. Therefore, the annihilation bound on 
particles with mx ^ MeV is only nonrcdundant for par- 
ticles that were not in thermal equilibrium at the time of 
nucleosynthesis . 
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Particles with mx ^ 5 MeV are only constrained by 
our bound if their density, extrapolated by (l + z)"'^ from 
the time they influence the CMB to the present, is sub- 
stantially greater than the present critical density, or if 
their cross section is very large but their number den- 
sity is somehow higher than the density obtained from 
freezeout of their annihilations (see Fig. |l|). The first 
case would require that the particle decays invisibly or 
otherwise disappears between z ~ 10^ and the present, 
and the second requires that the particle was formed by 
decays of a heavier particle (or some other mechanism) 
after the freezeout temperature for its annihilations. 

New particles that can be created in stars (e.g., by plas- 
mon decay, 7 XX) are constrained by their action as 
additional sources of cooling. Constraints from observa- 
tions of globular cluster stars apply for mx <10KeV|. 
Similarly, cooling in supernovae, specifically SN1987A, 
can be infiuenced by the creation of new particles with 
mass mx ^ 30 MeV; however, these constraints depend 
on the details of the particle interactions in the plasma 
(e.g., reference 39|). 

In summary, the bound from distortions of the CMB 
energy spectrum probably cannot be a useful constraint 
on annihilations of light {mx ^ 0.5 MeV) dark matter 
particles (e.g., warm dark matter) to photons, because 
the bound from annihilations in our Galaxy is always 
stronger. It is most interesting as a constraint on dark 
matter particles in the mass range 0.5 ^ mx ^ 5 MeV, 
and any particle that decays invisibly between z ~ 10^ 
and the present, although in either case the particle must 
also evade the bound from the expansion rate of the uni- 
verse during BBN. Finally, although we have not dis- 
cussed in this paper the case of annihilations of particles 
with cross sections that increase with decreasing temper- 
ature (e.g., [Q), it seems likely that they can be very 
tightly constrained by the kinds of tests we have dis- 
cussed. 
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FIG. 1. Bounds on relic particles from residual annihi- 
lations. The space above the lines is ruled out in each case. 
Except for the upper solid line, all bounds are for particles 
with {aa \v\) — (Jo (i.e., n — 0). The lower (upper) solid line is 
the constraint from the CMB (including both p and y distor- 
tions) for n = {n — 1). The upper (lower) long-dashed line 
is the constraint from BBN using the high (low) deuterium 
abundance (the low D constraint and the CMB constraint are 
practically identical but we have introduced a slight offset for 
clarity). The dashed line is the constraint from the diffuse 
photon background produced by extragalactic annihilations, 
and the dotted line is the constraint from the diffuse photon 
background produced by annihilations in the Milky Way halo. 
The BBN constraints apply for mx > 26 MeV, the thresh- 
old for photodissociation of ''He to D. For the CMB and BBN 
constraints, / is the fraction of the annihilation energy that in- 
teracts electromagnetically. For the photon background con- 
straints, / should be replaced by f-^, the fraction of annihila- 
tions that produce two photons with energy mx (the lines for 
these constraints are terminated at the electron mass because 
f-y « f is likely if other annihilation channels are open). 
For the case of equivalent particle and antiparticle, substitute 



^Ix, and CTo/6x 10 



ao/3xl0" 
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